Even-odd effects in finite Heisenberg spin chains 
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Magnetic superlattices and nanowires may be described as Heisenberg spin chains of finite length 
A'^, where is the number of magnetic units (films or atoms, respectively). We study antiferro- 
magnetically coupled spins which are also coupled to an external field H (superlattices) or to a 
ferromagnetic substrate (nanowires) . The model is analyzed through a two-dimensional map which 
allows fast and reliable numerical calculations. Both open and closed chains have different properties 
for even and odd A'^ (parity effect). Open chains with odd A'^ are known [S. Lounis et al, Phys. Rev. 
Lett. 101, 107204 (2008)] to have a ferrimagnetic state for small A'^ and a noncoUinear state for 
large N. In the present paper, the transition length Nc is found analytically. Finally, we show that 
closed chains arrange themselves in the uniform bulk spin-fiop state for even A'^ and in nonuniform 
states for odd A'^. 

PACS numbers: 75.75.+a, 75.30.Kz, 75.10.Hk, 05.45.-a 
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Antiferromagnetic (AF) chains of Heisenberg spins, 
when subjected to an external magnetic field (and possi- 
bly to a uniaxial anisotropy) , are known to arrange them- 
selves in a spin-flop (SF) state where neighboring spins 
are almost antiparallel and orthogonal to the field^ . How- 
ever, such a result is valid, strictly speaking, only in the 
thermodynamic limit. For a finite chain, boundary con- 
ditions and finite size effects are expected to induce modi- 
fications on the bulk spin-flop configuration, determining 
a non uniform canting along the chain. 

A one-dimensional (ID) classical planar model of a 
Heisenberg uniaxial antiferromagnet. 



T-Caf = ^ He cos{9i ~ 6*^+1) 

i=l 

N 

- '^{Hacos^ 9i + 2Hcos9^), (1) 



was introduced forty years agc^ to study a semi-infinite 
AF chain {N 00). In Eq. ([1]), He denotes the ex- 
change field. Ha the anisotropy field, and 9i is the angle 
that the magnetization of the z-th ferromagnetic layer 
forms with the direction of the external field, H (spins 
are assumed to be planar). The bulk SF phase appears 
for H > y/2HEHA + H\ {H > for zero anisotropy 
Ha). 

At that time, the reference experimental systems were 
bulk systems like MnF2 or MnO, with magnetic ions on 
special crystallographic planes interacting ferromagncti- 
cally (FM) and ions between planes interacting AF. With 
the spreading of epitaxially grown systems, the model 
was used a loli^i^i^i^iii^ to study superlattices made of N 
ferromagnetic layers which are antiferromagnetically cou- 
pled. Some important theoretical results were found: (i) 
For semi-infinite systems, the surface SF state (a phase 
predicted to anticipate the bulk SF state when increasing 
the field) does not exist'"^; (ii) For finite systems, there are 
important differences between structures with even and 
odd NM^. 



Recently, the model of a finite ID quantum Heisenberg 
antiferromagnet 
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has gained new interest since it has been used to describe 
an AF nanowire deposited on a thin insulating layer—. 
Paradigmatic examples of such a system are linear chains 
of 1 to 10 Mn atoms epitaxied on a CuN substrates'. From 
the analysis of spin excitations of coupled atomic spins 
in the dimer and in the trimer {N = 2,3), the Mn-Mn 
exchange interaction was found^ to be antiferromagnetic 
(|Ji| — 6.2 meV) and the spin value to be 5 = 5/2, 
identical to the spin of a free Mn atom. Using these 
parameters in Eq. ([2]), the magnetic behavior of longer 
wires could successfully be fitted^. 

When such AF Mn nanowires are deposited on a fer- 
romagnetic layer, like Ni(OOl), an interesting frustration 
phenomenon occurs, since the exchange coupling between 
an adsorbed Mn spin and the magnetic moment of an un- 
derlying Ni atom of the substrate is fcrromagnctioiSili, 
J2 > 0, and thus competes^ with the Mn-Mn antiferro- 
magnetic exchange, Ji < 0. Therefore, in a classical spin 
approximation, one is led to consider the modeUi 



n 



N 



N 

- J2 E '^OS 



(3) 



where 9i denotes the angle that the i-th spin of the AF 
nanowire forms with respect to the magnetization of the 
ferromagnetic substrate. Since the coupling J2 is local- 
ized on the i-th site, it is apparent that it plays the same 
role as the magnetic field 2H in Eq. ([1]), while | Ji| has to 
be identified with He, and Ha (the uniaxial anisotropy) 
is zero. Therefore, when N is finite and open boundary 
conditions are assumed, the existence of different ground 
states for odd versus even iV is a well-known result'^'". 
Different ground states also reflect on different behaviors 
for the spin wave excitations^. 
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FIG. 1: (color online) Phase portrait for mapping (|4)). Full 
red squares and blue circles correspond to equilibrium configu- 
rations for open chains with N = 9 and N = 10, respectively. 
Full arrows point to the values 8i for the first atom of the 
chains. Dashed arrows point to the values 9n+i (denoted by 
open symbols) for the fictitious (A*' -I- l)-th atom. Arrow tips 
allow one to follow the evolution of the map, whose hyperbolic 
and elliptic fixed points of order two, respectively denoted by 
SF and AF, correspond to uniform states. 



In a recent Letter— S. Lounis et al., using both ab 
initio results and solutions to the classical Heisenberg 
model ^ , confirmed that the ground state of finite AF 
nanowires deposited on ferromagnets depends on the 
parity of the number N of atoms. They also found 
that, while even chains always have a noncoUinear (NC) 
ground state, for odd N a transition from a coUinear fer- 
rimagnetic (FI) to a NC configuration occurs when the 
chain length N exceeds a critical value Nc. For example, 
using an iterative numerical scheme in order to minimize 
Eq. ([3]), the transition length was estimatedii to be 9 
atoms for Mn chains on Ni(OOl). 

Here we show that the classical Heisenberg model 
Q can be investigated with great numerical and ana- 
lytical profit in terms of a two-dimensional (2D) map 
method§''^>iii^ii^ii^ii^. Such an approach allows a fast 
and exact determination of the ground state configura- 
tion of finite chains and to find an analytical expression 
for the transition length for odd open chains. 

By the map method, we also study model ^ in 
the case of periodic boundary conditions. For such 
"closed" chains, we find a new even-odd effect: even 
chains have spins arranged in the spin-flop state, like 
infinite chains, while odd chains arrange themselves in 
noncoUinear states. 

In order to find the equilibrium configurations of the 
classical Heisenberg model ([3]) by the 2D map method, we 
introducei^ the variable s„ — sin(0„ — 9n-i)- Denoting 
by h = J2/IJ1I the ratio between competing exchange 
interactions, minimization of ^ gives^^ 

Sn+i = Sn - hsm0n, 0^+1 = On + sin"^(s„+i). (4) 

These equations define an iterative 2D map-S, i.e. 
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FIG. 2: (color online) Equilibrium configurations 6\,. . . ,6n 
calculated by the 2D map method for open chains with A'' = 
9, 10, 39, 40, and for a closed chain (PBC) with iV = 39. 



point [6, s) in the phase space is mapped to a point 
{6',s'). The fixed points of order two (s„+2 = Sn 
and 6n+2 = On) correspond to the coUinear AF con- 
figuration ((0,0) <-> (tTjO)) and to the bulk SF state 
((^,sin2^) ^ {-9,-am29), with cos^= h/4). In Fig.[T] 
we plot the fixed points and the evolution of the map for 
different initial conditions and h — 0.376 (it is the special 
value considered in Ref. [ll| as representative of AF Mn 
nanowires on Ni(OOl)). 

Boundary conditions for open chains of N atoms are 
taken into accounli^ii^ by introducing a fictitious (A + 1)- 
th atom and imposing si = = sat+i. The determina- 
tion of the ground state therefore corresponds to finding 
the value 6i such that, iterating the map N times from 
the point Pi = {Oi,0), we get a point Pn+i ~ {9n+i,0), 
with both Pi and Pn+i located on the horizontal axis, 
s = 0. The N values 6i, . . . ,6^ then give the sought-after 
equilibrium configuration. In Fig. [T] we also plot the first 
N steps of the map evolution giving the ground states 
for A = 9 (red solid squares) and = 10 (blue solid 
circles). Different behaviors for even and odd N can be 
inferred from the different location of their trajectories 
in the phase portrait. The configurations are explicitly 
shown in Fig. [2l 

The existence of a minimum length to get a non- 
coUinear configuration for odd N is clear from Fig. [31 
where we plot sn+i as a function of 9i, assuming si = 0. 
For odd A < 9, the only zeros are the AF fixed points, 
corresponding to a coUinear ferrimagnetic (FI) config- 
uration, but for A = 9, dsNJ^i/dOi changes sign at 
9i — 0, and an additional solution appears: the non- 
coUinear (NC) configuration. For even A, non trivial 
solutions exist already for A = 2 (dashed line). The In- 
set of Fig.[3]shows that for large N the function sn+i{9i) 
is strongly oscillating with several zeros 9^''^ . In order to 
determine the ground state, the energies of all the NC 
configurations with 9i ~ 9^'^'^ must be compared. 

We are now going to show that, by linearizing the map 
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FIG. 3: (color online) We plot sjv+i(^i), assuming si = 0, 
for different, small values of A'^ (main figure) and for A*' = 39 
(inset). The arrows point to the equilibrium values 9i for the 
first atom of the A'^ = 9 chain (red full line, full arrow) and of 
the N = 2 chain (dashed line, dashed arrow). 
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FIG. 4: Main: analytical results of the phase diagram for 
odd- A'' open chains, Eqs. (|11)I12|) . FI and NC denote Ferrl- 
magnetic and NonCoUinear states. The full line is the curve 
he — tt/N which is a very good approximation, even for small 
A' (Inset). 



nearby the fixed points, it is indeed possible to determine 
the exact analytical condition for the rising of the NC 
state in the case of an open chain with odd N. If we 
start from a point Pi = (si,0i) close to the fixed point 
(0, 0), iterated points on the map are osciUating between 
the AF fixed points: more precisely, even are close 
to (tt, 0) and odd P2k+i are close to (0, 0). So, if we write 

02k = T^ + hk (5) 

02k+l = 52k+l (6) 

the quantities J„ are small for any n, as s„ are. Now, we 
can linearize the map in the two cases P2fc-i — P2k and 
P2k P2k+i- If we write 



52k 
S2k 



Ax 



and 



S2fc+1 



A. 



hk 

S2k 



we get the matrices 



Ai = 



l + h -1 
-h 1 



and A2 = 



1 - h 
h 



-1 
1 



So, if 2iVo 

SN+1 



1 is an odd integer, we have 
\No ( \ = t( ^1 

Sl 



= AM2Aiy 



Since sm+i — 72i'5i + T22S1, if we start with si — 0, the 
condition dsN+i/dOi — dspf+i/d6i = reads T21 — 0. 
Let us now implement this condition, firstly determining 
eigenvalues and eigenvectors Vi of the matrix 



A = A2A1 = 
It is easily found that 



1 + h-2 
l-h 



Ai,2 = \{2-h'^ ±ih^/A-K^) 



Xj+h-l 
1 



1 



1,2. 



(7) 

(8) 



If U is the (2 x 2) matrix with vi^2 as column vectors 
and A]j is the diagonal matrix with elements Ai_2, it is 
straightforward to write T AiUA^"!/'^. Finally, the 
condition T21 = gives 



No 
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),No 
A2 , 



which simplifies to 



No 



A2- 1 

Ai-r 
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(10) 



Therefore, the transition length is equal to Nc — '2No + 
1 where Nq is the solution of the above equation. We get 



with 
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hV4: - 
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The curve is plotted as circles in Fig. [4] along with the 
asymptotic form Nc = Tr/h (full line) which appears to 
be a very good approximation even for small N (see the 
Inset). 

We now turn to closed chains, which imply periodic 
boundary conditions (PBC). If spins represent magnetic 
layers, these boundary conditions are not physical, but 
for nanowires deposited on a substrate they are phys- 
ical and correspond to nanorings. In terms of the 2D 
mapping, PBC imply Pn+i = Pi, i.e. On+i = Oi and 
sn+1 = Sl. Therefore, trajectories are fixed points of or- 
der N. It is easy to realize that the ground state for even 
N is the bulk spin-flop state: 92k — & and 92k+i = —9, 
with cos^ = /i/4. In fact, if 9i, . . . ,9^ were a different 
configuration with a lower energy, we might replicate it 
indefinitely for an infinite chain and get a configuration 
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FIG. 5: Graphical representation of the ground state for a 
closed chain with A'' = 9. Full and dashed arrows represent 
spin orientations on odd and even sites, respectively. 9i = 
and for every angle 92k there is an angle d2j+i = —82k- 

with an energy lower than the bulk spin-flop phase (which 
is the ground state). 

The above argument does not apply to odd N, because 
the SF phase, as well as the AF phase, which are fixed 
points of order two, do not satisfy PBC for odd N. In 
this case, we have a nonuniform state. In order to find 
it, we should look for the points P that are iterated on 
themselves after N applications of the map. It appears 
that the configuration is symmetric with respect to the 
field direction, i.e. for any spin with an angle 9* there is 



a spin forming an angle —9*, see Fig.[5]for = 9. There- 
fore, for odd N, there is one spin with 9 ~ Q. PBC allow 
to label this spin as "number 1" , so that searching the 
solution is now as easy as for open chains: we apply the 
map N times to points (0, si) and look for the values si 
such that 9n+i = and s^v+i — si- Using this method, 
we have found the ground state for = 39 (Fig. [21 full 
squares) and for = 9 (Fig. [S]) . 

In conclusion, we have studied model ([3]) which de- 
scribes a chain of classical planar spins with nearest 
neighboring AF coupling and interacting with a (real or 
effective) external field h. Since the anisotropy is zero, 
the infinite system is in the bulk SF phase for any nonva- 
nishing h. Extrema of the energy correspond to trajec- 
tories of the 2D mapping ^ with appropriate bound- 
ary conditions: si ~ sw+i = for open chains and 
^AT-i-i = 9i, S7V+1 — si for closed chains. 

This method allows a fast and exact determination of 
the ground states for any A^ (Fig. [2]). It also allows to 
find analytically the transition length Nc for open odd 
chains from the FI to the NC state (Fig. [5]). This tran- 
sition corresponds to a change in sign of the derivative 
dsN+i/d9i at 01 = (Fig. [3]). Parity effects are present 
for open and closed chains. 

With increasing the field h, the 2D map starts devel- 
oping a chaotic behaviori^^. Studying this regime with 
reference to nanowires would be an interesting subject 
for future work. 
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